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Abstract 

Observer design is considered for a class of non-linear systems whose non-linear 
part is energy preserving. A strategy to construct convergent observers for this class 
of non-linear system is presented. The approach has the advantage that it is possible, 
via convex programming, to prove whether the constructed observer converges, in 
contrast to several existing approaches to observer design for non-linear systems. 
Finally, the developed methods are applied to the Lorenz attractor and to a low 
order model for shear fluid flow. 

1 Introduction 

Observer design for non-linear systems is an important and difficult problem. In this 
paper, observer design is considered for systems whose non-linear part has an energy 
preserving structure. In particular, 

x{t) = Ax{t) + N{x{t))x{t), t > 0, 

y{t) = Cx{t), t>0, (1) 

x(0) = xo G M". 

where A £ M"^", C G M^^" and iV : M" — ^ M"^" is a linear operator. Furthermore, it is 
assumed that the non-linearity N(x)x has the energy preserving property 

x'^N{x)x = 0, X G M". (2) 

Observer design for non-linear systems has received much attention, with approaches 
falling into two main categories. One approach, first considered in and generalized 
in [HI [121 [111 ES]) is to apply a change of co-ordinates to linearize the system, up 
to an additional term involving the output y{t). Subsequently, linear design methods 
can be applied to create an observer for the transformed system, then the co-ordinate 
transformation is inverted to form an observer for the original, non-linear, system. The 
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main drawback of this approach is that it is usually impossible to prove that the chosen 
co-ordinate transformation is invertible. Hence, while this is a powerful technique for 
observer design, it is difficult to prove in practice that the constructed observer will 
actually converge. 

The second approach is to assume a Lipschitz-type bound on the non-linear part of the 
system. For example, in addition to the standard Lipschitz assumption [H [20], one- 
sided Lipschitz conditions [9l f23j and a 'less conservative' Lipschitz condition [19] have 
been studied. These techniques apply a Luenberger-type observer and require that the 
non-linearity is 'small enough' with respect to the linear part of the dynamics. A major 
drawback of this approach is that systems with a dominant non-linear term often have 
a large Lipschitz bound and, if this is the case, it is unlikely to be possible to prove that 
a given observer converges. 

The difficulties of the above techniques arise either from excessive generality or overly 
restrictive assumptions. The co-ordinate transformation technique may theoretically be 
applied to any non-linear system, and is therefore unlikely to succeed in every case. For 
the Lipschitz approaches, a small global Lipschitz bound restricts the class of systems to 
which the results may be applied. For this reason, we aim for an approach to observer 
design that sits between these two extremes by only considering the particular class of 
non- linear system ([T]) whose non- linear part satisfies Q. 

The importance of dynamical systems of the form Q-Q is that they often arise in finite- 
dimensional approximations of non-linear physical systems, for example, the Navier- 
Stokes equations in fluid flows [8j and the non-linear oscillation of beams in structural 
dynamics [7] . In experimental practice, such approximations are referred to as 'low order 
models' and can be created directly from experimental data by using, for example, the 
method of Proper Orthogonal Decomposition [10^ I21j. From a theoretical perspective, 
it is therefore of great interest to study the control theoretic properties of such systems, 
with a view to guiding experimental implementation. The link between the Navier-Stokes 
equations and ([l]) is presented in Section |3j 

Notation: The n-sphere in M"+-'^ is defined as S" := {x G M" | ||x||2 = !}• A matrix 
P G M"^" is said to be positive definite (written P ^ 0) if its symmetric part satisfies 
x^{P + P^x > 0, for any x G M", and negative definite if —P is positive definite (written 
P -< 0). The set of symmetric matrices of dimension n is denoted 5". For matrices A, B 
and C of appropriate sizes, the shorthand 



'A B' 
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P " 


. (*) c _ 




PT 


C 



is used to simplify the block matrix. For r > and d G M", the closed || • ||2-norm ball 
centered at d of radius r is denoted 



Br{d) :={xGM" I \\x-d\\2 <r] 



For i = 1, . . . ,n, 

Si = (0,...0,1,0,...,0)T 

^ V ' 

1 in i"^ entry 

denotes the i^^ element of the standard basis of M". For sets S,T C M" and a G M, 
S®T:={s + t:s£S,t£T}, aS := {as : s £ S}. 

2 Observer Design 

The approach taken in this paper is to exploit and energy preserving properties of the 
non-linearity in ([T]) to obtain a method for constructing a convergent observer. In par- 
ticular, for a given gain matrix L £ M"^P, the observer {x{t))t>o is assumed to have 
dynamics 

i = Ax + N{x)x-L{y-Cx), x(0) = xq G M". (3) 
Therefore, the observer error e := x — x satisfies 

e = {A + LC)e + N{x)x-N{x)x. (4) 

The aim of this paper is to find a constructive method of calculating L such that 

e{t) 0, t ^ oo. 



The main results. Theorem 2.6, Algorithms 2.11 and 2.13, and Theorem 2.15, provide 
methods of constructing such a gain L by solving a series of convex optimization prob- 
lems. 



2.1 A state invariant set 

The property of the non-linear system ([T]) that is advantageous for observer design is 
that the energy preserving property ^ implies the existence of an invariant set for the 
system dynamics. A set S C is said to be invariant for the dynamical system ([T]) 
if x{to) G S at time to > implies that x{t) E S for every subsequent time t > to- 
Invariant sets for the class of system ([T]) can be described in terms of perturbations of 
the linear part A of the system. In the following, given a matrix A G M"^" and a vector 
d G M" define a perturbed matrix Ad € M"^" by 

AdX := Ax + N{x)d + N{d)x, x € M". 

Subsequently, we make the following assumption. 
(Al) There exists d G M" such that Ad ^ 0. 



Clearly, assumption (Al) holds if ^ -< 0. Furthermore, it is shown in Lemma 3.1 that 
(Al) holds for the class of systems representing finite dimensional approximations of 
fluid flows. 



Lemma 2.1. Suppose that there exist d G M" and q > such that + al ^ Q. Then 
Br{d) is invariant for ([T]) for any 

r > -\\Ad + N{d)d\\2. (5) 
a 

Proof. Using the linearity of N : M" — t- M"^", the system's dynamics can be written 

x = Ax + N{x -d){x-d) + N{x)d + N{d)x - N{d)d 
= Ad{x -d) + N{x -d){x-d) + Add - N{d)d. 

For D{x) := ^\\x — d\\2, the energy preserving property ^ implies that 

D[x) = {x- dYAdix -d) + {x- dy{Ad + N{d)d). 
< -a\\x - d\\l + \\x - dhWAd + N{d)d\\2. 

Therefore, D{x) < whenever ||x — d\\2 > r, with r given by ([s]). Hence, B.,.{d) is 
invariant for ([T]). □ 

If ([T]) represents a fluid system, an invariant set may be calculated more explicitly, as 
described in Section 3.2 Since Ad is affine in d, the condition {d : + ^ 0} / can 
be checked by solving a semidefinite program [Ij. 

Ideally, one would like to calculate an invariant ball with the smallest possible radius. 
However, due to the non-linear dependence of (jsj) upon (a, d) G M x M", it is difficult to 
minimize ([s]) by convex optimization methods. In order to remove the non-linear depen- 
dence upon d from ([5]), the search can be restricted to vectors d such that N{d)d = 0. 

We first demonstrate that such vectors always exist. 

Lemma 2.2. Suppose that /:§"—)■ M"^^ is a continuous tangent vector field, that is 

x^f{x) =0, X G S", 
satisfying f{x) = f{—x). Then f has at least one zero on S". 

Proof. The result is well known when n is odd, in which case the condition f{x) = f{—x) 
is not required; see [22] for a particularly elegant elementary proof. 

We therefore consider the case where n is even. Assume that f[x) ^ for all x G S", 
and define the continuous function 5 : S" — )• S"^ as g{x) := /(x)/||/(x)||. Then g does not 
have a fixed point on S", since otherwise x~^ g{x) = x~^ x = 1 at the fixed point, which is 
not possible since g' is a tangent vector field by construction. 



Since g has no fixed points, its degree is odd ^ p. 134] (see [6, §2.2] for a definition of 
the degree of a continuous map). On the other hand, since g{x) = g{—x), its degree 
must be even [6, p. 165], a contradiction. Therefore / must have at least one zero. □ 

The existence of a nonzero x satisfying N{x)x = is then guaranteed by setting f{x) = 



N{x)x in Lemma 2.2, Since x i— >■ N{x)x is homogeneous, any such zero also satisfies 
N{ax){ax) = for all a G M. 

We henceforward make the following assumption: 

(A2) A matrix Q G M"^" is chosen such that N{d)d = whenever d G ker((5). 

The advantage of this assumption is that if the search for the centre of an invariant set 
is conducted over ker(Q), it can be performed by solving a semidefinite program. The 
existence of such a Q is guaranteed, since one can always define Q := I — zz^ where 



z 7^ is a zero of N{z)z whose existence is guaranteed by Lemma 2.2 



We note that identification of such a zero may be difficult in general. Define symmetric 
matrices Q^*^ such that Q^.^ := \\N{ej)ik + A^(efc)jj], for each i,j,k G {1, . . . ,n}. Then 



N{d)d= f(iTQ(i)d,...,(iTQ(")d 



(6) 



and computing a root of N{d)d amounts to finding a simultaneous root of n quadratic 
equations in n variables. See ^ i5j for numerical solution methods for such problems. 
However, if P|"^]^ ker(Q(*)) 13 {0} it is possible to select nontrivial Q such that 



ker(Q) = Pi ker {Q^'^ 



i=l 



We demonstrate the application of this method to the Lorenz attractor in Section [3j 
Even if this is not the natural choice for Q may be apparent given the system's 

underlying structure - see Section |3.2[ 



Proposition 2.3. Suppose that the semidefinite program 



subject to 



^ 



minimize s 

s {AzY 

(*) sin 

tA + {A^-A) + I^O 
t>0, Qz = 



(7a) 

(7b) 
(7c) 



with variables s,t £ M and z G M" has optimal solution (s* ,t* , z*) . Then t* > and 
Bs*{z* /t*) is an invariant set for {x{t))t>Q, with 

s* = inf {a~^\\Ad + N{d)d\\2 \ Ad + al ^0,d£ker{Q)] <oo. (8) 

a>0 



In the particular case s* = and z* = 0, then B^(0) is an invariant set for {x{t))t>Q for 
any 7 > 0. If ker An kerQ = {0}, then a minimizer to this SDP is guaranteed to exist 
if it is feasible. 

Proof. We first show that any feasible point of ([T]) has t > 0. Assume instead that there 
is some feasible point with t = 0, so that ( |7b[ ) satisfies {A^ — A) -< for some z. Then 
z^{Az — A)z = 2z^N{z)z < 0, which violates the energy conservation condition ([2]). 

We now consider the case when s* = 0. Let {t*,z*) be an optimal solution in this 



case. Then by rtTaJ, (7c) we have z* G ker A and N{z*)z* = Qz* = 0. Consequently 



(7b) guarantees that A + A~^ -< 0. In this case V{x) = \\x\\2 is a Lyapunov function for 



z* G ker(r^ + {Az* — A)) and (7b) implies that z* = 0. In particular A^* — ^ = and 



(^and a ball of any radius centered at the origin is invariant. 

If s* > 0, then ([T]) is equivalent to ^ after applying a change of variables a = 1/t, d = 
z/t and rewriting ( |7a[ ) as a quadratic constraint via Schur complement. The inequality 
( |7b[ ) is equivalent to + qI ^ (note the identity tA^/t = tA+ {Az — A)). Invariance 
then follows from Lemma |2.1[ 



The existence of a minimizer for ^ can be established by showing that the problem is 
equivalent to one with compact constraints. If ([7|) is feasible, then any value s = s at any 
feasible point can be used as an artificial upper bound on s. Since kervl n kerQ = {0}, 
there exists e > such that \\Az\\ > e||2:||, for any z G ker(^ Since (7a) is equivalent to 
||j42;|| < s, we can also add a constraint ||2;|| < s/e =: z. 

Denote as ^ the magnitude of the smallest negative eigenvalue of + A'^. Define 

C > as 

C := sup (Tmax / + l[{Az - A) + {Az - A)^ 

\\z\\<z I ^ 

Then one can also impose an upper bound t < t = C,/^ without altering the minimum 
value of ([T]). Augmenting the constraints in ([T]) with s < s, t < t and ||z|| < z, so that 
the feasible set is compact without altering theoptimal value, ensures the existence of a 
minimizer. 



Remark 2.4. Note that the kernel constraint (7c) is included in order to ensure that 
the problem ([T]) is solvable as a semidefinite program, and is conservative in the sense 
that it restricts the search for an invariant set Br{d) to those with centers satisfying 
N{d)d = 0. 



^ Suppose for a contradiction that there exists a sequence z„ € kerQ \ {0} such that < 
^11 2,1 II, n > 1. Let Zn/\\zn\\- Then 5„ G and hence there exists a convergent subsequence 

^n,. -)■ C G S""^ n kerQ. However, 

\\A^\\ = lim PCnJI < lim — = 



which imphes that ^ e ker j4, contradicting ker A n kerQ = {0}. 



□ . 



It is also possible to remove this condition and solve the more general problem 



inf \a-'^\\Ad + N(d)d\\2 



(9) 



directly, i.e. to solve the problem ^ without a kernel constraint. Assuming that s* > 0, 
one can make a change of variables a = 1/t, d = z/t and apply the Schur complement 
to get the equivalent problem 



minimize 



subject to 



s {Az + N{d)zy 

(*) sin 

tA + {A^- A)+I <{) 
t > 0, z = td. 



>- 



(10a) 

(10b) 
(10c) 



Noting that z = td is equivalent to the pair of constraints diag(z) < diag(td) ^ diag(2;), 



the constraints in optimization problem ( 10 ) constitute a set of bilinear matrix inequali- 



ties (BMIs). Although methods for solving optimization problems of this type are avail- 
able ]13^ . there is generally no guarantee that a solution will be globally optimal. We 
therefore follow the somewhat more conservative method of Proposition \2.3[ 



2.2 Locally stable observers 



If it is possible to calculate an invariant set for the state, Theorem 2.6 of this section 
provides a strategy for constructing a locally convergent observer. First, it will be useful 
to derive an explicit expression for the norm of the non-linear term A^. 

Lemma 2.5. Suppose that N is given by ([g]). Define matrices Q*-^^ = {qlj^)^j^i by 

\\N\\2 = CTmax(6)2, 

where Q is a matrix whose {i,jY^ entry is @ij := {q'^^\ Q^^'^^ , and {A,B) = tr{A^B) 
is the standard (Frobenius) inner product. 



Proof. It can be shown from ^ that N{x) = YA=iXiQ^^ , for x = (xi)^^i G W. The 
adjoint N* : M"^" W is given by 



)W\ 



iV*(Q) = QeM"^", 

and hence, N*N = Q. □ 

Theorem 2.6 (Local Observer Convergence). Suppose that Br{d) is invariant for {x{t))t>o 
and y C M"" satisfies Br{d) C Y. Suppose that there exist {ai)f^i > 0, P £ and 
R £ ]R"^P such that 

aiI<P<a2l (11) 



and 

PAy + + RC + C^R^ ^ -asl, 
If {x{t))t>o has dynamics ^ for L := P^^R and 

as lai 



|eo||2 



\xq - X0II2 < 



(12) 



(13) 



2.5 



27Q2 V ai 

then e{t) — )• 0, t — )• 00. The constant 7 := (Tmax(0)^ is defined in Lemma 
Proof. Using 

N{x)x — N{x)x = N(x — x)x — N{x)x + N(x)x 
= N(x — x)x + N{x)(x — x) 
= N{e)x + N{x)e - N{e)e, 

the error dynamics Q can be rewritten 

e = {A^^t)+ LC)e- N{e)e, t>0. 

Since the underlying state dynamics {x{t))t>o are unaffected by (e(t))(>o it is possible 
to consider the time varying linear operator A{t) := A^^^^ + LC independently of the 
error dynamics. Hence, if V{e) := e^Pe, 



(14) 



V(e) = (PAit) + A{typ) e - 2e~^PN{e)e 



(by ([TT]), ([T2J)) < -a3||e||2 + 2a2||e||i||A^(e)||2 
(by Lemma [23]) < (-Q3 + 27a2||e||2)||e||2. 



(15) 



Now let e G (0, 03) be such that ||eo||2 
(11) and (15) imply that, 

V{e) < aip^ =^ I 



Q3 — € / Ql 



2702 V 012 



and define p := (03 — e) 72702. Then 



||e||2 <P^V{e) < -e\\e\\l 
By (11) and (13), V{eo) < a2||eo||2 = aip^. Hence, {V{e{t)))t>o is decreasing and 



< 



Vieit)) ^ F(eo)„-^ 



< 



e < pg 



0, 



t — 7- 



□ 



Remark 2.7. A simple, but instructive, necessary condition for {11), (12) to hold is 
that the pair (A^, C) is detectable. In other words, the output map must at least be 
compatible with the linear system generated by the perturbed matrix A^. Furthermore, if 
By{d) is invariant for (x(t))(>o with A^ + otl^Q and r = a^^'^Ad + N{d)d\\2, then 



< --d^ Add = --d'^ {Ad + N{d)d) < -I 
a a a 



\Ad + N{d)d\\2 = r\ 



Hence, G Br{d). Therefore, in the case that the state invariant set is calculated by 
Proposition^^ detectahility of {A,C) is also necessary for (11), (12) to hold. 

A sufficient condition for local convergence can be formulated involving only the matrices 
Ad and C . If there exists a > and ^ P ^ a{A'yr)~^ such that 

PAd + ^JP + RC + C^R^ ^ -al, 



then it is not difficult to show that (12) holds for Y = Br{d), implying that the observer 
is locally convergent. 



We discuss a method for reformulating the semi- infinite LMI constraint (12) as a finite- 
dimensional LMI in Section (2.5). 



2.3 An observer invariant set 



If a locally convergent observer can be constructed by Theorem 2.6, it is natural to ask 
whether it is possible to extend the set of initial states for which the observer converges. 
Since it is known that the system state {x{t))t>o has an invariant set S, say, it is desirable 
for the observer {x{t))t>o to itself possess an invariant set which contains S. 

The following two results provide a method for calculating an invariant set for the 
observer dynamics. The first of these characterizes the trapping set for the observer 
error dynamics, and parallels the results of Lemma |2.1[ 



Lemma 2.8. Suppose that Bj.{d) is invariant for ([T]) and that xq G Br{d). Suppose that 
there exist d £ and a > such that A^ + LC + al ^ 0, then Bf{d) is invariant for 
(|3| for any 

f>- sup LCv - Ad- N{d)d . (16) 

'-^ V&Br{d-d) ^ 

Proof. Suppose that d £ W is such that A^ + LC + al <0. Then, 

i = {A^ + LC){x-d) + N{x - d){x - d) 

- {LC{x - d) - Ad - N{d)d). 

For D{x) := ^\\x — d\\2, the energy preserving property ([2]) implies that 

D{x) = {x- ciy{A^ + LC){x - d) 



(x - dy{LC{x -d)-Ad- N{d)d) 



< -a\\x-df 



2 

{x - dy{LC{x - d) - Ad - N{d)d). 



By assumption, x{t) G Br{d), for any time t > 0, so that 

x{t) - dG Brid-d), t>0. 



Therefore D{x) < whenever \\x — d\\ > r, with f given by (16). Hence, Bf.{d) is 
invariant for ([s]). □ 

The next result provides a method for computing an invariant set for the observer dy- 
namics given an observer gain L. As in case for the state invariant set, the non-hnear 
dependence of (16) upon (a, d) G M x M" makes global minimization difficult. We there- 



fore remove the non-linear dependence upon d from (16) by restricting the search to 



vectors d such that N{d)d = 0. In the following, the assumption is made that the ob- 



server gain L is such that A-\-LC is stable. Note that, by Remark 2.7, this is a necessary 



condition even for local convergence. The following result parallels Proposition 2.3 



Proposition 2.9. Suppose that Br{d) is an invariant set for ([T]) with r > 0, A + LC 
stable, and that the semidefinite program 



minimize 



subject to 



s-rt\\LC\\2 {LC{td- z) - AzY 
(*) (s-rt||LC||2)/„ 
t{A + LC) + {A^- A) + I <^ 
t > 0, = 



^ 



(17a) 

(17b) 
(17c) 



with variables s,t G M and z £ M" has optimal solution (s* ,t* , z*). Then (s*,t*) > 
and Bs*{z* /t*) is an invariant set for {x{t))t>o, with 



s* > inf <^a-^sup \\LCv - Ad - N{d)d\\2 A^ + LC + al ^ 0,d ekei{Q) \ . 



(18) 



A minimizer to this SDP is guaranteed to exist if it is feasible. 



Proof. If LC = 0, then A is stable and the result follows from Proposition 2.3 



Now suppose that LC 0. Using the same argument as in the proof of Proposition 2.3 



one can show that any feasible point of (17) must satisfy t > 0. The constraint (17a) 
then requires s > rt\\LC\\2 > 0. 

Noting that Qd = =^ N{d)d = by assumption, the smallest invariant set radius f 



satisfying the inequality (16) for a given a simplifies to 
f= sup a~"^ ||LCi; — Aci||2 

veBr{d-d) 

= sup a~'^\\LCv + LC{d- d) - Ad\\2 

VGBr{0) 

< ra'^\\LC\\2 + a-^WLCid - d) - Ad\\2- 



(19) 



Upper bounding (19) by s and substituting t = 1/a and z = d/a, results in 

f < rt\\LC\\2 + \\LCtd - (A + LC)z\\2 < s. 
and applying a Schur complement identity produces the equivalent linear matrix inequal- 



ity (17a). One may likewise confirm that the inequality vl; + LC + a/ ^ is equivalent 



to (17b). Invariance then follows from Lemma 2.8 



To establish the existence of a minimizer, we first show that the SDP (17) is equivalent 
to one with compact constraints. Any value s = s at any feasible point can be used as 
an upper bound on s, which allows an additional constraint t < t := s/ {r\\LC\\2) to be 



imposed as a necessary condition for (17a). A further necessary condition for (17a) is 
then 

||(^ + LC)z||2 < s + \\LCd\\t. 

Noting that kei{A + LC) = {0} since (^4 + LC) is assumed stable, the remainder of the 
proof proceeds as in the proof of Proposition |2.3[ 

□ 



Remark 2.10. Note that the proofs of Propositions 2.3 and 2.9 are similar, hut that the 



result of Proposition 2.3 produces a tight hound on the invariant set radius for the state 



dynamics ([T]), whereas the result of Proposition 2.9 is conservative due to the application 
of the triangle inequality in (19). 



2.4 Globally stable observers 



To study observer convergence, it is useful to rewrite the nonlinear part of the observer 
error dynamics (|4]). Starting from (14), for x, a; G M" the nonlinear error term is 



N{x)x - N{x)x = N{e)x + N{x)e - N{e)e 

= N{e){x - e/2) + N{x - e/2)e 

= N{e){{x + x)/2) + N{{x + x)/2){e) 

= {A^ - A)e. 



Hence, the error dynamics (|4]) can be written 



Ax+i + LC) e. 



(20) 



(21) 



The observer error dynamics can therefore be considered as a linear time varying system. 



and the problem of observer design is to find a gain L which stabilizes (21). Since both 



the state and error dynamics can be contained inside separate invariant sets Br{d) and 



Bf{d) respectively, our objective is to identify a gain L that stabilizes (21) under the 
assumption that 

{x{t) + x{t)) eBr{d)®Bf{d), t>0. (22) 



The central difficulty is of course that the estimation error trapping set Bf[d) is itself 
determined by the observer gain L. We therefore propose a two-phase strategy, which 



we characterize formally in Algorithm 2.11 



Our general approach is first to identify a trapping set Bj.{d) for the state dynamics ([T]) 
using the method of Proposition 2.3 We then select some set Y such that Br{d) C Y, 
and compute a gain L such that Ay + LC is stable for all y G Y. Using this gain, one 
can compute an observer invariant set Bf.{d) using the results of Proposition 

0. 



2.9 



If such 



a set exists and {Br{d) © Bf{d))/2 C y, then (21) is stable and e{t) 
Algorithm 2.11 (Observer design). 

1. Use Proposition 2.3 to select d G M",r > such that Br{d) is invariant for ([T]). 



2. Select (01,02) > and Y C M" such that Br{d) C Y. Compute a positive definite 
Pes'" and Re R^'^'p such that: 



PAy + AlP + RC 



P-ail y 

a2ln RC 
(*) a2lp 

C^R" 



0; 
0; 

< 0, 



yeY. 



(23) 
(24) 
(25) 



Define L := P-^R. 



3. Use Proposition 2.9 to select d G M",f > such that Bf{d) is invariant for ([s]). 

4. If \ Br{d) ® Bf{d) C Y, then ([21]) is stable and e{t) — )■ 0, whenever xq G Bf{d). 



Remark 2.12. The tuning parameters (01,02) appearing in Step^ of Algorithm \2. ll\ 
are included to provide control over ||LC||2. Recalling (19) in the proof of Proposition 



2.9. \\LC\\2 influences the radius of the observer invariant set calculated in Step^ Min- 



imizing the size of the this set is useful in helping to ensure that the set inclusion in Step 
[2 is satisfied. 

To ensure that Proposition \2.9\ can be used to construct an invariant set for the observer 
in Algorithm 2.11 , Step^ one must of course first verify that 



jo o > 0,3(i G M",^j + LC + a/ -< 0} / 



The fact that R and P are searched for simultaneously in (23)-(25) may mean that 



IILCII2 is suboptimal, and consequently that Step [4] of Algorithm 2.11 does not hold. 
In this situation, we propose the following iterative search for a globally convergent 
observer. 



Algorithm 2.13. Initialization : Suppose that Steps\l^ -[^ of Algorithm \2. ll have been 
completed to provide: 

1. d G M", r > such that By{d) is invariant for ([T]); 



2. PeS'',R€ M"^P and Y C M" such that Step\^ of Algorithm \2Al\ holds; 

3. d £ M", r > such that Bf{d) is invariant for Q. 

Define Pq := P, Lq := P^^R, do := d, f-Q := f, oq := and (3o := 0. 



Iteration : repeat until ^ Br{d) © Br^{dk, 



C Y: 



1. Suppose that a* > and P* E 5" are an optimal solution to the semidefinite 
program 

maximize a 
subject to P y al 

P{Ay + LkC) + {Ay + LkCyP yeY. 

Let Pk+i := P*,ak+i := a* . 

2. Suppose that /3* > and L* G M"^^ are an optimal solution to the semidefinite 
program 



minimize /3 



subject to 



>- 



/3/„ LC 
(*) pip 

Pk+i{Ay + LC) + [Ay + LC)^Pk+i < 0, 



Let Lk+i :=L*,l3k+i :=/3*. 

with L 
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3. Apply Proposition 

Bf^^-^{dk+i) is invariant for ([s]). 



Lk+i to select dk+i G W^,fk+i > such that 



4. If \ Br{d) © Bf^^-^{dk+i) C Y, then e{t) for any xq G (dfc+i)- 

Remark 2.14. At each stage of the iterative proceedure, {Pk, ^miniPk)) is feasible for 
the SDP in Step^, while (L^, HLfcClb) is feasible for the SDP in Step^ Furthermore, 
for each k >1, we have the bound \\LkC\\2 < /3k/ak- 



An alternative to the iterative method of Algorithm 2.13| is to search for P over a 
particular subset of 5", defined in terms of the non-linearity N. If this subset is well 
chosen, it is possible to remove the need to find Bf.[d). Define 



{P G 5" : e''PN{e)e = 0, for each e G M"}. 



Since the energy preserving property Q holds, it is the case that 6^ 7^ 0. Notice also 
that, since e^PN(e)e is linear in P, it is easy to calculate SJ^ for a given non-linearity 
N. The following result provides conditions for global observer convergence. 



Theorem 2.15. Suppose that By.{d), calculated by Proposition 2.3, is invariant for ([T]) 
and let xq € Br{d). Pick y C such that Br{d) C Y and suppose that there exists a 
positive definite P E SJq and R € M"^^' such that 



PAy + AyP + RC + C'R' -< 0, 



yeY. 



Then if L := P ^R, the observer {x{t))t>o defined by ^ satisfies e{t) 
any initial condition xq € M". 



Proof. Let V{e) = e^Pe. Then by (20) and (21) 



VV-e = e'^P{A^ + LC)e 

2 

= e^P{A + LC)e 

+ e''P{N{x)e + N{e)x - N{e)e) 
{P G = e^P{A^, + LC)e. 
(by ^) < 0, e G M". 



(26) 

0, t — )• oo, for 



Hence, e{t) — 0, i — >• (X) and since x does not appear in the expression for W ■ e, the 
observer error converges to zero for any initial condition xq £ M". □ 



2.5 Modeling of Robust LMI Conditions 



In order to apply the results of Theorem 2.6, Algorithms 2.11 and 2.13, or Theorem 



2.15, it is necessary to construct matrices P and R such that the semi-infinite matrix 
inequality 

PAy + aIp + RC + C'^R^ ^ -al, yeV, (27) 

is satisfied for some compact set Y C and for some a > 0. We next comment on 
methods for modeling such a constraint as a finite-dimensional LMI. 



Suppose that Y = Br{d), so that (27) can be rewritten as 

PAd + AjP + RC + C^R^ + PiAs -A) + {As - A^ P ^ -al, 
Define ^ := {P, R,a), a matrix F(o)(0 G 5" as 



< r. 



FW(0 : = 



PA, + A\P + RC + C^R^ +al 



) 



and matrices F^^\(^) G for i G {1, . . . ,n} such that 

n 

-P{As -A)- {As - A)'P =: ^ J,F«(e). 



i=l 



The robust LMI condition (27) can be rewritten in this notation as 

n 

F{0)(^) + ^5^i.W(^)^0, ||<5||2<r, 



i=l 



where each of the matrices F^'^\^) is hnear in ^. We can then exploit a result from robust 



semidefinite programming to establish a sufficient condition for satisfaction of (27). 



Proposition 2.16 ([31 Thm. 2.1]). IfY = Br{d), then the robust LMI ([27j) is satisfied 
if there exists a > 0, Q G 5" and 5 G 5" such that {S + Q) ^ 2F(°)(0 and 



y 0. 



In the more general case that Y C conv{?/j,i = 1,...,M}, one can of course also 



guarantee satisfaction of the constraint (27) by ensuring its satisfaction at every vertex 



3 Examples 

We give two examples of observer design for finite dimensional systems related to fluid 
flows; the Lorenz attractor and a low order model for shear flow between two parallel 
plates. 



3.1 Lorenz Attractor 

The Lorenz attractor [T7] is a dynamical system in M'^, which is a simplified model of 
fluid convection in two spatial dimensions. We consider the classical Lorenz dynamics 
which can be written in the form ([T]) for 



A 



N{x) := -xi , X G M^. 
\0 XI / 

Note that non-linearity N satisfies the energy preserving property Suppose that it 
is possible to observe only the second state: 



-10 


10 


28 


-1 





- 











-Xl 


XI 







y = Cx, C := ( 1 ) 



A convergent observer can be constructed using the iterative method of Algorithm 2.13 



Proposition 2.3 implies that Br{d) is invarient for {x{t))t>o with 



r := 100.7, 



d 



37.5 



Select a conservative boun ding set Y = Brg{do) with ro := 1200 and do := d. After five 
iterations, Algorithm 2.13 provides matrice^ 



-10.0 



-13.3 



Ps = diag( 1000.0 -0.1 -0.1 



which satisfy P^iAy + L^C) + {Ay + L^C)^ ■< 0,y £ Y. The set Bf^{d^) is invariant 
for the observer dynamics ([s]) with 

h = 1282.6, = ( 9.2 



and it can be easily verified that ^ Br{d) ® Bf,^{d^) C Y. Hence, Algorithm 
implies that x{t) — )■ x{t),t — )• oo whenever 

{xo,xo) G Br{d) X Bf^Xd^). 
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Alternatively, Theorem 2.15 can be used to construct a globally convergent observer. 
For the Lorenz attractor, we have 

5]^ = spanjdiag ( 1 ) , diag ( 1 l)}. 



and it is interesting to note that the matrix P5 constructed by Algorithm 2.13 



is an 



element of S^. Applying Algorithm 2.11 , Step[2] with the restrictions P E 5]^, P -< 10^1 
provides matrices 



L = ( -9.6 -704.4 



P = diag( 132.4 0.8 0. 



€5" 



N 



(28) 



2.15 



implies 
i^. An 



which satisfy P{Ay + LC) + {Ay + LC) ^ P 0,y € Br{d). Hence, Theorem 
that the resulting observer is globally convergent for any initial value xq G 
example of the performance of the two globally convergent observers is shown in Figure 

m 



Observer design for the Lorenz attractor is considered in [16], where the co-ordinate 
transformation approach is used. This approach creates an observer which appears to 
converge experimentally, but the complexity of the co-ordinate transformation means 
that it is not possible to prove convergence. 



For the Lipschitz approach, suppose there exists 7 > 0, P G 5" and R G M*'^" such that 

-< 0. (29) 



PA + A^P + RC + C^E7 
(*) 



P 



■^The additional condition P < 10''/ was imposed to improve convergence. 




It is easy to deduce (see e.g. [9l [19]) that if S is an invariant set for the state and the 
non-hnearity satisfies the Lipschitz condition 

||iV(x)x - N{y)y\\2 < 7||x - y\\2, x,y £ S, 

then ([3]), for L = p-^R, is a convergent observer. With respect to the Lorenz dynamics, 
the largest 7 > satisfying (29) is 7 = 2.67. However, letting x = (xi, X2, X3), y = 
{yi,X2,X3) implies that 

\\Nix)x-N{y)y\\2 _ 



\x - y\\2 



X2 I /Y"2 
2 ' '''S- 



It is known that there exists x in the range of the Lorenz attractor for which y^x|~+~r| > 
\/l500 and hence, the Lipschitz approach (e.g. from [T|) cannot be used to construct a 
convergent observer for the Lorenz attractor. 



3.2 Low order model for shear fluid flow 

We consider observer design for the finite dimensional fluid flow model presented in 
|18j . The model is derived from the Navier-Stokes equations by the method of Galerkin 
projection. Before considering the example, we explain how this method necessarily 
results is a system of the form ([T]) with nonlinear term satisfying the energy preserving 
property Q. 



The incompressible Navier-Stokes equations for a vector field u : Q C 



, are 



1 

— + {u-V)u = -Vp + -^V^u + /, 
ot Re 

V-u = 0, 

where p : — )• M represents the pressure, / : 17 — )• M'^ an external force and Re the 
Reynold's number of the flow. No-slip boundary conditions {u^g^ = 0) are also assumed. 
A common assumption |10^ [2T] is that the flow field can be decomposed in the form 



u{x,t) = ai{t)ui{a 



(30) 



i=l 



and a finite dimensional approximation of the flow obtained by considering the truncation 

En 
i^i aiUi. 

A set of ordinary differential equations for the time-dependent coefficients Oj can be 
obtained via the method of Galerkin projection (see e.g. fSj pp. 129-154]), leading to 



{f,Ui) Xi 



Re 



1, . . . ,n. 



(31) 



where Aj > are fixed constants. To remove the constant term from (31) it is assumed 
that there exists a known stationary point a = c. Making the transformation x = a — c, 
the perturbations about c have dynamics of the form ([T]) with linear part 



Ax 



Re 



Ax + N{c)x + N{x)c, X e M", 



(32) 



for A 



-diag ( Ai • • • A„ ) ^ and nonlinear term 

n 



i=l 



As a consequence of the incompressibility and no-slip assumptions, 

{{uj ■ V)uk, Ui) = -{{ui ■ V)uk, Uj), 

which implies that the matrices Q^*^ are anti-symmetric. Hence, the nonlinearity N 
satisfies ([2]). 

As a first step towards designing an observer, we construct an invariant set for the state 
dynamics. Although Proposition 2^ can be applied, the particular structure of the linear 
term (32) implies that a natural invarient set can be easily constructed. 



Lemma 3.1. Suppose that {x{t))t>o satisfies (1) with linear part of the form (32), for 
some c € M", and nonlinear part satisfying (2). Then there exists r > such that 
Br{—c) is invariant for {x{t))t>o- 



Proof. Note that A_c = -< 0. By Lemma 2.1 it follows that ||x + cHl is decreasing 
if 

-^(x + c)^A(x + c)T+ {x + cY ( -^Ac + N{c)c] < 0. 
Re \ Re J 

Standard algebraic manipulation shows that the set of x E for which the above 
inequality holds is equal to \ E, where E is the ellipsoid 

2 



£; := |x G M" : A, (^Xi - ^ (Ac + 7V(c)c),^ 

i=l * 

The result follows if r > is chosen such that E C Br{—c). □ 

We now consider observer design for a low order model of shear fluid flow. For brevity, 
we refer to |18l pp. 7-8] for an explicit description of the modej^ and note that all 
subsequent calculations are performed for Reynolds number Re = 60. 



For this system, the first vector field ui appearing the expansion (30) coincides with 
the laminar solution to the flow, implying that c = ei. Since A^(ei)ei = 0, Lemma 3.1 
implies that ei) is invariant for the system if 

In particular, Bx^/\^^^{—ei) = Bi{—ei) is an invariant set. In fact, applying Proposition 
2.3 with kerQ := span(ei) implies that 

^^(-Cei), i := 0.9477 



^ := <( X G M" : Ai ( xi + - ) + > , A^x,^ < ^ K Br{-ei) 



is an invariant set. Hence, Proposition |2.3| provides a tighter invariant set that than the 
natural one derived from Lemma |3.1[ 



We assume that the first six states of the system can be observed, i.e. 



where O G M*^^^ has all entries equal to zero. Let y be a 1-norm ball of radius 3.^ such 



that B^{—^ei) C Y . Applying Algorithm 2.13 with 10 iterations provides an observer 
gain Lio and Pio G 5" for which 

Pw{Ay + LioC) + {Ay + LioC)Pio ^0, yeY. 



Consequently, Theorem |2.6| implies that the observer is locally convergent. The com- 
plexity of the system makes it unlikely that a globally stable observer can be constructed 
using our methods. However, Figure [2] shows the unobserved states of the locally stable 
observer can be seen to converge to the true system state. 



'With respect to the system parameters in [18], we select a — 1/2, /3 = 7r/2 and 7 = 1. 



3;(i) f(i) 




t t 



Figure 2: The state fluid system with initial condition 

(-0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5)"^ is denoted {x{t))t>o. The observer with 
gain LiQ is {x{t))t>o with initial condition (0 0000000 0)^. 

4 Conclusions 

A method of observer design has been presented for a class of non-linear systems whose 
non-linearity is energy preserving. Sufficient conditions, which can be verified by stan- 
dard convex optimization methods, are given which imply either local or global observer 
convergence. The results are applied to create a globally convergent observer for the 
Lorenz attractor and a locally stable observer for a low order model of shear fluid flow. 
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